Astronomy & Astrophysics manuscript no. 1116ms 


February 2, 2008 


(DOI: will be inserted by hand later) 





Electric field screening by a proton counterflow in the pulsar 

polar cap 



Katsuaki Asano 1 and Fumio Takahara 2 



^1" 

o 
o 

(N 
< 



> 
On 
00 

o 

O 

^ : 

i ' 
o ■ 

a: 

C/3 . 



1 Division of Theoretical Astronomy, National Astronomical Observatory of Japan, 2-21-1 Osawa Mitaka Tokyo, 
Japan 

e-mail: asano@th.nao.ac.jp 

2 Department of Earth and Space Science, Osaka University, Toyonaka 560-0043, Japan 
e-mail: takahara@vega. ess . sci . osaka-u. ac . jp 

the date of receipt and acceptance should be inserted later 

Abstract. We propose a new mechanism to screen the electric field in the pulsar polar cap. Previous studies 
have shown that if an electron beam from the stellar surface is accelerated to energies high enough to create 
electron-positron pairs, the required electric field parallel to the magnetic field lines is too strong to be screened 
out by the produced pairs. We argue here that if non-relativistic protons are supplied from the magnetosphere 
to flow towards the stellar surface, they can provide an anode to screen out such a strong electric field. Injected 
electron-positron pairs yield an asymmetry of the electrostatic potential around the screening point. The required 
pair creation rate in this model is consistent with the conventional models. 
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1. INTRODUCTION 

A spinning magnetized neutron star provides huge 
electric potential differences between different parts 

of its surface as a result of unipolar induction 

ijGoldreich fc Julian Ul96^1 . A part of the potential differ- 
ence may be expended on an electric field along the mag- 
netic field somewhere in the magnetosphere. Although a 
fully self-consistent model for the pulsar magnetosphere 
has not yet been constructed, several promising mod- 
els ha ve been proposed. Among them, the polar cap 
model l|Sturrock L Il97lt iRuderman fc Sutherland! Il975|) 
assumes that an electric field E\\ parallel to the magnetic 
field lines exists just above the magnetic poles. The elec- 
tric field accelerates charged particles up to TeV ener- 
gies, and resultant curvature radiation from these par- 
ticles produces copious electron-positron pairs through 
magnetic pair production. These pairs may be sources for 
gamma-ray emission, coherent radio emission, and the pul- 
sar wind. 

The localized potential drop is maintained by a pair 
of anode and cathode regions. In the cathode region 
the space charge density p deviates from the Goldreich- 
Julian (GJ) density negatively. On the other hand, p 
must deviate positively for the anode. Outside the ac- 
celerator the electric field will be screened. In the po- 
lar cap model, especia l ly for a space charge-limited flow 
model l|Fawlev et al. L Il977t IScharlemann et al. L Il978l 



lArons k. Scharlemann I Il979|) , where electrons can freely 
escape from the stellar surface, i.e., En — on the stel- 
lar surface, the formation mechanism of a static anode- 
cathode pair that can sustain enough potential drop for 
pair production is a long-standing issue. Since the current 
flows steadily along the magnetic field lines, the charge 
density is determined by the magnitude of the current 
and the field geometry with suitable boundary conditions. 
Basic considerations on t his problem for space charge- 
limited flows are given in IShibata ] Jl997h . The mecha- 
nism of the electric field screening, i.e., a way to provide 
an anode, may be provided by pair polarization. Although 
most papers take it for granted that copious pair produc- 
tion can instantly screen the field, recently Shibata et al. 
(1998, 2002: hereafter SMT98 and SMT02, respectively) 
have casted doubt on this issue; electric field screening is 
not as easy as most researchers consider conventionally. 

SMT98 and SMT02 investigated the screening of elec- 
tric fields in the pair production region. They found that 
the required thickness of the screening layer is restricted 
to be as small as the braking distance Ze — TO cC 2 /l e -^j| L 
for which decelerating particles become non-relativistic, 
where m c and — e are the electron mass and charge, 
respectively. If the above condition does not hold, too 
many positrons are reflected back and destroy the neg- 
ative charge region (cathode). In order to screen the elec- 
tric field consistently, a huge number of pairs should be 
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injected within a region as small as 1~e- The required pair 
multiplication factor per primary particle is enormously 
large and cannot be realized in the conventional pair cre- 
ation models. Thus, some other ingredients are required 
for the electric field screening. 



In relation to this problem, iLvubarskii I 1^)92) pro- 
posed a model with a positron countcrflow and argued 
that the model makes it possible to screen the electric 
field for any current density. He showed that the electric 
field on the stellar surface can be screened by the injection 
of electron-positron pairs plus a positron counterflow from 
the magnctosphere for any current density. In his model 
the positron counterflow is introduced to adjust the cur- 
rent desnity in the polar cap and does not play a critical 
role as an anode formation. Most importantly, his model 
is different from the space charge-limited flow considered 
here and cannot be applied to our problem. As SMT98 
and SMT02 have shown, the way to screen the electric 
field at a reasonable altitude is unknown for space charge- 
limited flows which we believe to be the most suitable 
model because the work function of the matter on the 
neutron surface is considered to be small enough. 

In previous studies of the screening, pairs were as- 
sumed to accelerate or decelerate along the 0-th order 
trajectories determined by E\\ . However, if an electro- 
static (longitudinal) instability occurs, the excited waves 
may produce effective fricti on which will change the 
charge polarization process. lAsano k. Takaharal <j2004^l 
numerically investigated an electrostatic instability of 
electron-positron pairs injected in an external electric 
field. However, they found that the density of electron- 
positron pairs in the standard pulsar model is too low to 
induce the instability needed to provide strong friction. 
This means that each particle moves stably along the tra- 
jectory determined by the external electric field since col- 
lective forces from other particles are negligible. 

In this paper, we present a new model with a proton 
counterflow to solve the screening problem in the polar cap 
model with pair injection. Since the proton mass is much 
larger than the electron mass, its inertia effect is much 
stronger and can provide an anode structure more eas- 
ily if non-relativistic protons are supplied from the outer 
magnctosphere to flow towards the polar caps. Our pur- 
pose is to point out a possibility that the proton flow can 
explain the structure of the electric field in the polar cap. 
Therefore, our model is highly idealized and simplified. 
Detailed investigations of the situation in our model are 
beyond the scope of the present paper. In Sect. 2 we briefly 
summarize the space charge-limited flow model in the pul- 
sar polar cap to give a basic framework for the present 
model. In Sect. 3 we consider the model with proton coun- 
terflow and demonstrate the effectiveness of the screening 
by the proton counterflow. Sect. 4 is devoted to a summary 
and discussion. 



2. THE POTENTIAL FOR SPACE CHARGE 
LIMITED FLOW 

Before we present the model with a proton counterflow, 
we summarize the basic features of the behaviour of the 
electric potential in the pulsar polar cap to make clear 
the fundamental points of the problem. The notations are 
the same as in SMT98 and SMT02 in which the one- 
dimensional (ID) approximation is used. We assume that 
SI* B > 0, where O* and B are the angular velocity of the 
star and magnetic field, respectively. Although the current 
density distribution is most likely to be determined by the 
global dynamics in the magnetosphere, we will investigate 
a local region just above the stellar pole. We therefore 
treat the current density as an adjustable free parame- 
ter in our paper, as assumed in SMT98 and SMT02. A 
steady electron beam flows from the stellar surface and is 
accelerated by a negative electric field parallel to a mag- 
netic field line with zero electric field on the stellar sur- 
face. Beyond a certain altitude called the pair production 
front (PPF), electron-positron pairs are assumed to be in- 
jected continuously into space. In this case a part of the in- 
jected positrons returns towards the stellar surface. Here, 
we consider mainly the behavior of the electric potential 
inside the PPF while that outside the PPF was treated in 
SMT02. 

Since the current due to the primary electron beam 
flowing outward and that due to positrons returning to- 
wards the stellar surface flow inward along the magnetic 
field, both currents are negative. Inside the PPF the dif- 
ference in the charge density determined by these currents 
and the GJ density gives the effective charge density. We 
define the GJ current density as 

We normalize the potential length along the magnetic 
field line s, GJ current density Jo, current density due to 
the primary electron beam J, and current density due to 
the returning positrons J+ as 



Jo 



m c c 

LOpS 
C 

Jb' 
J 



J+ 

3+ = -j-, 

J B 

where 



Jb 



2^r 

Aire | Jb 



<0, 



(2) 
(3) 
(4) 
(5) 
(6) 

(7) 

(8) 



The current densities J, J+, Jo, and Jb are all negative, 
so that jo, j, j+ are all positive, and by definition jo < 1. 
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The equation of continuity ensures that j oc J/B and 
j+ oc J+/B are constant along the magnetic field line. On 
the other hand, jo °^ E z /B, where B z is the component 
of the magnetic field along the rotation axis, may change 
along the field line. When the magnetic field line curves 
away from the rotation axis, jo decreases with increasing I. 
Conversely, when the field line curves toward the rotation 
axis, jo increases with I. Adopting the ID approximation, 
the non-dimensional Poisson equation becomes 



d 2 ^ 



L + 2± 



Jo 



J_ 

3 



J+ - Jo, 



(9) 



where /3 = v/c > and (3 + — v +l c < are velocities of the 
primary electron beam and the returning positrons nor- 
malized by the light velocity, respectively. Since positrons 
are ultra-relativistic inside the PPF, we have approxi- 
mated as (3 + ~ — 1 here. 

The boundary conditions are given on the stellar sur- 
face (I = 0) as follows; 4> — 0, d(f>/dl = 0, and the electron 
beam is supplied with (3 — 0. In this case, the first term 
due to the electron beam, on the right hand side of Eq. 
(jHJl, dominates the other terms on the stellar surface. The 
non-relativistic electron beam provides a huge negative 
charge density at I — 0. As a result, just above the stellar 
surface d 2 fbjdl 2 > 0, which means that the potential curve 
is downward convex, irrespective of j+ and jo- As long as 
the first term j / (3 dominates, <b increases with increasing 
I and the electric field becomes negative. The electric field 
accelerates the electron beam, and (3 becomes ~ 1 at a 
certain altitude. Hereafter, we represent the scale of the 
transition from (3 ~ to ~ 1 as l nr ; for I > Z nr , we write 
3 ~ 1. For I > l nl , the Poisson equation is approximated 
as 



dl 2 ~ 
where 



-j+j+ +3o, 



(10) 



(11) 



is the effective charge density inside the PPF except for 
< I < l nr . In general, the scale l m is much smaller than 
the length from the stellar surface to the PPF. In the 
greater part inside the PPF, the potential obeys Eq. I|10|) . 

When we assume that jo is constant neglecting effects 
of field line curvature, p- m is constant. In this case, the 
behaviour of the potential is determined only by pi n and 
the situation is classified into the two cases: p- m < and 
Pin > 0. Hereafter we call the case of p; n < (j > j + + jo) 
Super-GJ and that of p m > ( j < j + + jo) Sub-GJ. 

2.1. Sub-GJ 

First, we discuss Sub-GJ. In the region inside Z nr , the non- 
dimensional electric field E = —d<f)/dl < grows with I 
as was mentioned before. For Sub-GJ, d 2 (j>/dl 2 for I > l nr 
is negative, so that the potential curve is upward convex. 
In this region the absolute value of E declines with I as 




Fig. 1. Potential for Sub-GJ and constant jo- The position 
I = i nr is the inflection point. For I < l nr , the effective 
charge density is negative, owing to the non-relativistic 
electron beam. For I > l nr , the effective charge density 
becomes p m > 0. 

schematically shown in Fig. 1. The electric field can be 
screened at a certain distance even without pair formation. 
In the approximation of constant p in , the solution of Eq. 
HI L)| l is obtained as 



(12) 
(13) 



E — — — —E + Pi n l, 

dl 

where Eq > represents the magnitude of the electric 
field at I = 0. Strictly speaking, for I < l nl the above 
solution is not correct and the electric field grows from 
at / = to — Eq inside l nl . The length l nr is negligibly small 
when considering the region / > l nl . From Eq. I|13[l . the 
screening altitude, where the electric field becomes zero 
again, is approximated as 

U = ^. (14) 

Pin 

If there is no pair injection, the electric field is not 
screened but changes sign at l scr so that the electrons are 
decelerated outside / scr . Thus, the model becomes incon- 
sistent. Let us further assume that around Z scr electron- 
positron pairs are injected and the electric field is screened. 
We represent the Lorentz factor of the beam at / = l scr 
as r(= 0(Z scr ) + 1) = 10 6 r 6 . Using T, the electric field at 
I = is written by 



E ~ V2p~T. 



(15) 



Substituting this into Eq. 1)14)1 . the altitude of the screen- 
ing is obtained as 



(16) 
(17) 




5 x 10 4 r 1 /V 1 / 2 T 1/2 /^ 1/2 

O X 1U 1 6 p in oJ-o.3 D X1 
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where p ln _ 3 = p m /lCT 3 , B 12 = B/10 12 G and T = 
0.3T .3 s is the rotation period of pulsar. 

To estimate Z scr , we must estimate plausible ranges of 
Eo, pi n and T. The value of T is considered to be 10 6 ~ 10 7 
to obtain copious pair production. To constrain Eq and 
Pi n , let us consider the behaviour of the potential just 
above the stellar surface (I < l m ). Multiplying Eq. © by 
E and integrating it from the stellar surface to an arbi- 
trary height, 



,1=1 



(18) 



where we use the relation Edl = —dcj>. The above condition 
corresponds to the Gauss law; the column charge density 
between I and the stellar surface is equal to E(l). Then 
we obtain 



- \e 2 = Oo + 3+)<t> - jV0 2 + 20. 

For Sub-GJ, \E\ takes the maximum value at 
which is given by 



(19) 



0m + 2 <j>n 



.4 2 



I - A 2 ' 



where A = j / (jo + j+) < 1. From Eq. 1)190. we obtain the 
maximum electric field as 



1 



,Ei 



= (1 - y/1- A 2 )(j Q + j+) < JO + j+, 



(21) 



We consider that this electric field gives an estimate of Eq 
in Eq. I|15|) . Then we have 



El 



2ph 



< 



jo + 3+ 

Pin 



(22) 



Since T is very large, pi n should be very small as long as 
jo and j+ are of order of unity. 

In order to complete the electric field screening, the 
effective charge density should be zero outside Z scr which 
is identified with the PPF. Outside the PPF, excess pair 
electrons provide the normalized current j- which is equal 
to the current j+ due to returning positrons inside the 
PPF. Using the approximation that both electrons and 
positrons are relativistic, the effective charge density out- 
side the PPF is given by 



Pout 



-J - 3+ + 3o = 0. 



(23) 



Thus, j+ = jo — j should be satisfied for self-consistency. 
The value j+ to be positive requires j < jo , which is equiv- 
alent to the Sub-GJ condition without pair injection. Then 
the constraint of a very small p- m = 2(j Q — j) = 2j+ im- 
plies that the number of returning positrons is very small 
and that E is of order of unity (see Eq. iJHJ). If we as- 
sume r ~ 10 6 , pi„ ~ 1CP 6 , which implies a high screening 
height of s scr ~ 10 6 cm from Eq. (TTJ. This altitude is 
much higher than conventionally supposed and actually, 
geometrical effects will play a role either in screening or 
in changing the flow into a Super-GJ flow. 




Fig. 2. Potential for Super-GJ and constant jo- The effec- 
tive charge density is negative everywhere independently 
of the velocity of the electron beam. 



(20) 2.2. Super-GJ 



For Super-GJ, since pi n < and d 2 (f>/dl 2 > the potential 
curve is downward convex everywhere (see Fig. 2), The 
potential and the amplitude of the electric field steadily 
increase with I. The electric field is given by Eq. I|13(l and 
becomes 



E ~ p in l, 



(24) 



for large I by neglecting Eo- If we conventionally assume 
that pair production can screen the electric field and if 
we also assume that the screening height l SCI roughly cor- 
responds to the PPF, the screening height turns out to 
be given by Eq. I|17|). when p m is replaced by |pi n |- The 
non-dimensional electric field around the PPF is obtained 



as 



En, 



V2|pin|r = -44.7 x 



Pin, 



,|l/2 r 
-3 1 J-6 



1/2 



(25) 



If | /3; n | is large, E scr becomes large around the PPF, 
while the screening height becomes low. As was shown 
in SMT02, pair polarization cannot screen the field when 
-Escr is more than a few. If E scr is lower than a few, the sit- 
uation is the same as for Sub-GJ, i.e., geometrical effects 
play a role without pair production. Another problem of 
Super-GJ is the following. For Super-GJ, the condition 
that the effective charge density outside the PPF vanishes 
becomes 



Pout 



-J - 3+ + 3o = Pin - 2j+ = 0, 



(26) 



where all the components are assumed to be relativistic. 
This condition is never satisfied for positive j+ within the 
ID approximation. In order to achieve screening for Super- 
GJ, some components should be non-relativistic. The ex- 
istence of the anode implies that the effective charge is 
positive at the point where E = 0. Therefore the effective 
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charge density should decrease outside this point. As a re- 
sult, the models may require that a space charge density 
wave continues outside the PPF as was shown in SMT98. 
In this wave, the speeds of electrons and positrons go and 
return repeatedly between a non-relativistic speed and the 
light speed. The average effective charge density outside 
of the PPF is zero, while the local charge density oscillates 
spatially. 

3. THE PROTON COUNTERFLOW MODEL 

Our proton counterflow model is motivated by the diffi- 
culties discussed in the previous section and provides a 
mechanism for realizing relevant E SCT and i scr for a rea- 
sonable T. As was shown in SMT98 and SMT02 and dis- 
cussed in the previous section, it is difficult to screen the 
electric field in the system composed of the electron beam 
and injected pairs only. For Super-GJ, we need to screen 
an electric field that is as strong as given by Eq. I)25|l. In 
this section we show that non-relativistic protons flowing 
backward can screen the electric field owing to the large 
proton mass. 

3.1. Origin of the Proton Counterflow 

To resolve the screening problem, we add a proton coun- 
terflow as a minimal alteration to the standard model. The 
existence and origin of the proton counterflow is highly 
speculative and it is difficult to give definite mechanisms 
at the present stage. Although it is highly unlikely that 
protons flow towards the polar cap from outside the light 
cylinder, the proton counterflow can be realized if pro- 
tons exist in the magnetosphere. Considering that we still 
do not well understand the fundamental issues of pulsar 
physics such as the current closure problem, we should be 
free from any kind of prejudices. 

There are some candidates for the origin of the pro- 
ton counterflow. One of the most plausible candidates for 
the source is protons existing in the co-rotating magneto- 
sphere. In the 0-th approximation such protons are con- 
fined to closed magnetic lines so that the protons do not 
flow into the polar cap region. However, if transport phe- 
nomena like the anomalous diffusion found in laboratory 
plasmas occur, the protons can flow into the polar cap 
region across magnetic field lines. Electrostatic (longitu- 
dinal) drift-waves are considere d to be the mo st dominant 
diffusion process in tokamaks l)Liewer I . llQSfJl . while vast 
amounts of theoretical efforts have been devoted to the 
understanding of the anomalous diffusion. 

A density gradient of particles is likely to exist on the 
boundary between open and closed field lines. The drift 
motion due to the density gradient occurs along the di- 
rection (hereafter y-axis) perpendicular to both field lines 
and the density gradient (x-axis) . Electrostatic drift- waves 
are excited accompanying the drift motion, so that the in- 
duced electric field along the y-axis induces E x B drift 
and particles diffuse along the x-axis, which will provide 
the source of the proton counterflow. The diffusion coef- 



ficient for this process may be obtained from numerical 
simulations. 

The source of the anomalous diffusion is not limited 
to the electrostatic drift-waves due to the density gradi- 
ent. The extremely high brightness temperature of pulsar 
radio emission requires a coherent source for this radia- 
tion and excitation of strong plasma waves in the open 
field line region. These waves may affect the corotating 
magnetosphere to induce proton diffusion. 

As an example, let us roughly estimate the amount of 
protons that diffuse from the co-rotating magnetosphere. 
Assuming that the diffusion occurs around the "outer re- 
gion" (radius Ro ~ -Rl/2, where Rl = c/f2* is the radius 
of the light cylinder). Since the magnetic field is relatively 
weak (typically I0 3 -I0 6 G), the drift velocity v± ~ cE/B 
due to the excited electric field can be large enough. Here, 
we assume that the scale of the proton diffusion region AR 
is much smaller than Ro- We denote AR = fRo (/ <C 1) 
hereafter. Then, the number of protons escaping from the 
co-rotating magnetosphere is written as 



2TrR ARSn v 1 _ = 7r/(c/0*) 2 <5n u±/2, 



(27) 



where Sno is the proton number density resonating with 
the excited waves. 

As will be shown below, our model requires that the 
proton current density on the polar cap is of the order 
of the GJ-current density. The typical size of the polar 
cap region i? po i on the stellar surface is ~ 10 4 cm. In our 
model the number of protons flowing into the polar cap 
region per unit time is 



2^ ' 



(28) 



where is the magnetic field on the stellar pole. The 
supply rate given by Eq. I|27|) is larger than the required 
rate if 



on > 2 2 , — 
it c ej v± 

~ 700 x 



T 0.3 B 12 



( R 



pol 



Vl0 4 c 



-Yi—Y 1 

0.1 VO.lc/ 



(29) 



cm" 3 . (30) 



This value seems very small for the atmosphere around 
such a high density celestial object, and it is less than the 
Goldreich- Julian density at Ro, though there are large 
uncertainties for v \ and /. In some atmosphere models 
l|Ho fc Lai I. l200ft IZane et al. I l200l[) that are consistent 
with observed proton-cyclotron resonance in a soft gamma 
repeater or an anomal ous X-ray pulsar ijlbrahim et al. I . 
2003: iRea et "Zl l2003h . the proton density is assumed to 
be larger than 10 20 cm~ 3 . Although the situation here is 
quite different from such models, the proton density in 
the magnetosphere may be much higher than the value of 
Eq. H30f) . Therefore, we consider that it is not so unlikely 
that a fraction of the protons diffuse from the co-rotating 
magnetosphere flow into the polar cap region, and that the 
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Fig. 3. Schematic picture for the proton counterflow 
model. 



global dynamics in the magnetosphere adjusts the current 
density distribution to screen the electric field. 

Further detailed discussion of the diffusion depends on 
the wave excitation mechanism, the model of global cur- 
rent flows, the model of the atmosphere, and the convec- 
tion in the co-rotating magnetosphere etc., which are not 
well understood. We will not discuss those problems in 
this paper. There may be mechanisms for the diffusion 
from the co-rotating magnetosphere other than described 
here and we do not stick to a specific mechanism. The 
existence of protons may be connect with the current clo- 
sure problem and the source and role of proton flows are 
interesting subjects in connection with all the problems in 
pulsar physics. 



We assume that the PPF is located in region 1 or 
II. Positrons in regions I and II are decelerated, and a 
part of the positrons returns towards the stellar surface. 
Therefore, in most of region I there are three relativistic 
flows: the electron beam coming from the stellar surface, 
the returning positron flow, and the proton counterflow. 
In this region the effective charge density determined from 
these flows is assumed to be negative (Super-GJ). As is 
discussed in Sect. 2.2, the amplitudes of the potential and 
electric field grow with I in region I. As the potential in- 
creases with /, protons become non- relativistic and their 
speed decreases with I. As a result, the positive charge 
density due to protons increases. In regions II and III the 
effective charge density is positive and the potential curve 
is upwards convex, because of the contribution of the non- 
relativistic protons. In region IV the effective charge den- 
sity once again becomes negative, because the speed of 
protons increases with I. 

The bulge of the potential is formed by the positive 
charge density due to the proton counterflow. Electron- 
positron pairs cause an asymmetry between region I-II and 
region III-IV. If there is no pair injection, region I-II and 
region III-IV are symmetric and the valley becomes <fr = 0. 
In this case, the potential becomes periodic within the ID 
approximation. When pairs are injected, as is required for 
pulsar action, part of the positrons cannot cross the po- 
tential peak and return towards the stellar surface. On the 
other hand, all injected electrons cross the peak and flow 
into region III. Thus, there is an asymmetry of the charge 
density between regions I-II and III-IV. In our model the 
local effective charge density in region IV (~ is much 
larger than that in region I (|pi n | "C 1). The existence of 
excess negative charge due to pair electrons in region IV 
leads to an immediate change of d<f>/dl with I, and will 
keep the valley shallow enough. 



3.2. Overview of the Model 



3.3. Region of E < 



The schematic picture for the proton counterflow model 
is depicted in Fig. 3. As shown in this figure, our model 
incorporates a bulge of the potential, which is divided into 
the following four regions: region I (E < 0, effective charge 
density p cS < 0), II (E < 0, p eS > 0), III (E > 0, p cS > 
0), and IV (E > 0, p e ff < 0). The potential increases 
from the stellar surface through regions I and II and peaks 
at the boundary between regions II and III. At the peak 
of the potential the electric field becomes zero. Then the 
potential decreases and the electric field becomes positive 
in regions III and IV. Our model requires that the electric 
field vanishes again in region IV. The potential curve in 
region IV is downwards convex and reaches the bottom 
again (we call this point "valley", hereafter). In our model 
we assume that pair electrons do not return, so that the 
potential drop should be smaller than the electron energy 
at injection. As was mentioned in Sect. 2.2, the potential 
continues to oscillate outside of the valley. The potential 
curve shown in Fig. 3 is possible if the proton counterflow 
exists, as will be explained below. 



In this subsection we discuss regions I and II. Following 
SMT02, outside the PPF the potential satisfies the 
Poisson equation, 

d 2 4> 



dl 2 J p p 



Jo 



P-, 



(31) 



where j p > and j3 p < are the normalized current 
density and velocity of the proton counterflow, respec- 
tively, and p + = cp + /\Jb\ and p~- = cp-/\Js\ are non- 
dimensional charge densities of pair positrons and elec- 
trons, respectively. Since pairs are injected for a finite 
length; their densities are calculated by integrating over 
the injection place, which will be done below. 

Let us set cj> = and E = E F at the PPF. At the 
peak of the potential (0 = 4>m), namely the boundary 
between regions II and III, the electric field becomes zero. 
Integrating the Poisson equation between the PPF and 
the peak as was done in Sect. 2.1, we obtain 

d<P(j - jo) - / d4>(p+ + P-) 



2 F 



/ 

Jo 
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„ d4> W 



Jp 



(32) 



The first term of the right hand side in Eq. (|32[1 becomes 



#(i -jo) = (j - Jo) 



(33) 



To calculate the second term we follow the same pro- 
cedure as SMT02. We introduce the multiplication factor 
m(l) in order to take into account the contribution of in- 
jected electron-positron pairs. The net pair flux produced 
between the PPF and a position I is denoted by m(l)j. 
Since </> and m are both monotonic functions of I in this 
region, m or <fi can be used to designate the coordinate 
position. The charge densities of pairs injected by a small 
flux element between m and m + dm are given by 



dm 

dp± = ±J W\ = ± 



J7± 



71-1 



zdm 1 



(34) 



where the Lorentz factors 7± are measured at the posi- 
tion (p, and are functions of the injection point denoted by 
^> in j(m), i.e., 



7_ = cj) - inj (m) + 7 
for electrons and 

7+ = -</>+ 0injM +7 



(35) 



(36) 



for positrons; 7 is the Lorentz factor at injection, and as- 
sumed to be constant for simplicity. 

A detailed analytic calculation of the second term in 
Eq. H32fl in general cases is shown in SMT02. Here, for 
simplicity, we adopt a toy model for the pair injection 
rate; the pair creation rate is constant with respect to 
the '^-coordinate', as dm = a^dcj), where is constant. 
While oicf, is the pair multiplication factor created in a unit 
potential drop m e c 2 /e, taking a constant is a good ap- 
proximation for a spatially uniform pair injection model 
(see SMT02). In this case the charge density of pair elec- 
trons is easily calculated as 



P-{4>) 



-3<*<t> 



(37) 



At a given position 0, the Lorentz factor of positrons 
injected at 4>- m j = — 7 + 1 becomes 7+ = 1. These 
positrons turn around there and return toward the stel- 
lar surface. For <f> > 7 — 1, positrons injected between 



7 + 1 < 



< 



move forward at a position 



Positrons injected at 



< 1 



7 return before 



<fi = 0m and do not enter region III, while positrons in- 
jected at 0i n j > 1 + 4>m ~ 7 enter region III and do not 
return. Therefore, at a given position <f>, positrons moving 
backward were injected at (/> — 7 + 1 < 4>- m - s < 1 + (j>M — 7- 



Thus, the charge density of pair positrons at <f> is given by 

1 



P+WO 



J 



, + i J 1/3+1 

1+0M-7 
+jOt < j > / d0inj 
-7+1 



1 



\0, 



(38) 



where the first term of the right hand side is due to the 
positrons moving forward, and the second is due to the 
positrons moving backward. For <fi < 7 — 1, we should 
change all the lower bounds of the integrals to zero. As a 
result, we obtain 



P+ = ja<t> S 



V7 2 " 1 + V(^m 
(V^^+V^M-^+l) 2 -! 



+ 1) 2 - I 

for <t> > 7 — 1, 



(39) 



-2V(7-0) 2 -l 
for (j) < 7 



1. 



The charge density of pairs is straightforwardly inte- 
grated to give the second term in Eq. (|3"2|) as 



d4>{p+ + pJ) = ja^ (2</> 



■b 

1 

-2 ( 



+ i)V(<I>m +7) 2 - 1 



(PM i- 

-- log 

2 B 



1)VW>M + 1) 2 -1 

G(0m+7)G(7) 



G(0 M + 1) 



where 



G{4>) = 4> + \f¥^- 



(40) 



(41) 



The third term of Eq. I|32(l represents the contribution 
of protons. The Lorentz factor of protons is given by 



7p = {4>m - (t>)8 + 7m, 
where 6 = m c /m p and 7 n 



(42) 

1 is the minimum Lorentz 



factor of the proton counterflow at <j) = 4>m- The third 
term in Eq. I|32(l becomes 



d ^Wl = -f(V(Hu+lm) 2 1 - !)• (43) 

\Pp\ ° 

Since we have calculated the righthand side of Eq. (|32(l , 
we now investigate whether the proton counterflow can 
make E-p large enough. One constraint on the parameters 
is imposed by the flow inside the PPF, because we are 
treating Super-GJ with p- m < 0. Since all components are 
relativistic inside the PPF except very near the stellar 
surface, the charge density is approximated as 



jo -j+j++j P < 0, 



(44) 



where j + is the non-dimensional current density due 
to returning positrons inside the PPF. Inside the PPF, 
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positrons injected at < 0i n j < 1 + 0m — 7 are returning. 
The charge density of positrons inside the PPF (0 < 0) is 
obtained as 

P+ = 3<H ( V(<Pm - 4> + 1) 2 - 1 - V(i - <P) 2 - 1) • (45) 

By definition, we have the relation 

3+ =ja<^(l + 0M ~7)- (46) 

It is seen that inside the PPF the charge density of 
positrons is approximated by p + = j + with high accu- 
racy. 

First of all, let us confirm that Ep should be smaller 
than the value in Eq. I|25ll when there is no proton coun- 
terflow. In this case 

r<f>M 

El/2= / d^(p+ + p-) - (j - j^u. (47) 
Jo 

In the limit of 0m > 7 > 1 and A = 0m - 7 > 1, we 
obtain 

El/2 ~ j+(0 M /A + 7) - (j - io)0 M 

= (j--j o )(0m/A- A)+p in (0 M /A + 7). (48) 

We should notice that j — jo > and p ln < for Super- 
GJ considered here. The second term on the righthand 
side is always negative. If A ~ 0m (namely 0m 3> 7), the 
first term ~ — A(j — jo) becomes negative too. Therefore, 
positive Ep is possible only for A 2 <C 0m- As is understood 
from Eq. |15|l. Ep cannot be larger than ~ 10 3 , for A > 1, 
3 - jo ~ O(l) and M ~ 7 = 10 2 -10 3 . If A is of the order 
of unity or less, Ep can be large. In this case, however, 
a0 = j+/j(l + A) becomes of the order of unity, while 
CV0 <C 1 in realistic situations (see SMT02). Thus, only 
Ep of the order of unity can be screened through pair 
polarization. 

When a proton counterflow exists, the situation 
changes drastically. The third term of Eq. I|32|l is easily 
calculated, and we obtain 

El/2~ ^ v /(,50M+7m) 2 -l+J+(0M/A + 7) 


-(j-jo)0M, (49) 

where the first and second terms are due to the proton 
counterflow and pairs, respectively. As long as the proton 
flow is non-relativistic, 7 m is not important for determin- 
ing Ep, though the local charge density depends strongly 
on 7 m . In conventional pair creation models, 7 is assumed 
to be a few hundreds, namely 7 ~ 0(1/ S). If 0m 3> 1/5, 
El/2 ~ (p- m - j+)(f>M + 3+1 = /3in0M - j'+A < 0. Too 
large a value of 0m leads to a large negative contribu- 
tion of the last term. Therefore, Ep becomes positive and 
large only when 0m, 1/5 and 7 are of the same order, 
but 1/5 is a few times larger than the others; in this case 
the contribution of the second and third terms becomes 
~ — Aj-f. — J P 0M- Even if A ~ 0m, this contribution is as 
small as ~ — O(j p 0m) for j+ ~ j p . The proton contribu- 
tion is a few times j p /6 and can provide a large value of 



Ep ~ y/jp/6. When j p is of the order of unity, Ep is a 
few times ten, which is close to the value in Eq. (|25|l . The 
small value of d makes Ep large enough. When j + 3> j p , 
the contribution of — Aj+ in the second and third terms 
dominates the other terms. Therefore, the current j + is at 
most comparable to j p . 

3.4. Region of E > 

In this subsection we discuss regions III and IV, where 
the electric field becomes positive so that positrons are 
accelerated while electrons are decelerated. The speed of 
incoming protons increases with I. There are two possibili- 
ties to complete the screening; one is to vanish charge den- 
sity and the other is to obtain a small-amplitude charge 
density wave with vanishing average charge density, which 
was shown in SMT98. The first possibility is unlikely be- 
cause this requires that the effective charge should vanish 
suddenly at the point where E — 0, although the charge 
density (proton density) tends to increase toward the peak 
(p e fi > 0) of the potential and decrease toward the valley 
(pes < 0). Thus, the second possibility is more natural. 
The oscillation amplitude should be small enough since if 
the amplitude in the potential is too large in this region, 
injected electrons would return to region II and continue 
to oscillate around the potential peak without escaping to 
infinity. 

In this subsection we discuss the condition for produc- 
ing the valley. The Poisson equation was given in Eq. Q31|) . 
Part of the pairs injected in regions I and II flow into re- 
gions III and IV. The injection region for such electrons 
is < < 0m in regions I and II, and for positrons it 
is 1 + 0m — 7 < < 0m- As for pairs injected in regions 
III and IV, we assume dm — ~a^d4>, because decreases 
with / in regions III and IV, Then we obtain 

P- = 3<*<t> ( 2 V (7 + <t> - M 2 - 1 - V1 2 - 1 

-V(7 + *) a -l), (50) 
p+ = ja (2V(7-^ + M 2 - 1 - V7 2 "! 

-V / (^m-0+1) 2 -1). (51) 

Since our model requires that the oscillation amplitude of 
the potential is small enough, all electrons and positrons 
injected in regions III and IV are relativistic. Therefore, 
the charge densities of such electrons and positrons cancel 
out, and scarcely contribute to the total charge density at 
all. The charge density due to pairs (p_ + p+) becomes 
~ — in regions III and IV, which is confirmed from the 
above equations. 

In region III the positive charge due to protons domi- 
nates over the negative charge density due to pairs, be- 
cause of the low speed of the proton counterflow. As 
the speed of the proton flows increases with I, the neg- 
ative charge due to pairs dominates in region IV and 
the valley will appear. Further out the spatial wave of 
the potential and charge density will repeat; this will 
be discussed in the next subsection. The condition to 



Asano and Takahara: Electric Field Screening in the Pulsar Polar Cap 



9 



make the effective charge density negative in region IV 
is ip/IA>l < 3+ + U - Jo) = 2j+ + j P - Pin- On the other 
hand, opposite inequality applies in region III because of 
the positive effective charge density in region III. A small 
value of /3 p in region III satisfies the condition. 

The condition in region IV is rewritten as \f3 p \ > 
j p /(2j+ + j p — pi„). In order to obtain a larger value of 
Ef, a smaller value of j+ <^ j p is favorable, while a fi- 
nite value of j+ is required to create the assymetry in 
the potential. From the above conditions, (3 P should vary 
significantly (~ 0(1)), which requires that the typical am- 
plitude of potential variation is of the order of 1/8. At the 
same time, pair electrons and positrons are assumed to 
flow outward continuously, which requires a fairly large 
value of 7 >^ 1/8 to overcome the potential variation. 

In the next subsection we will confirm numerically 
whether the model can reproduce both a shallow valley 
and a strong Ef- 

3.5. Numerical Demonstration 

We here numerically demonstrate the screening by the 
proton counterflow model in order to confirm the be- 
haviour of the potential in regions III and IV especially. 
In this subsection we assume that j = 1.0, p m = —0.001, 
7 = 500, and /3 m = -y/l- 1/7^ = -0.1. The other pa- 
rameters are j p , j + , and the potential difference between 
the PPF and the peak, </>m- The physical quantities j and 
are determined by the above parameters. We assume 
that the currents j p , j+, and j are all 0(1). From Eq. (|32|l 
we can determine the electric field at the PPF, Ef- By 
setting E = Ef and <j> = at the PPF, we can numer- 
ically solve the Poisson equation from the charge densi- 
ties we have obtained in this section. In Table 1 we show 
some examples of our numerical results that succeed in 
screening the electric field with subsequent charge density 
wave. As discussed in Sect. 3.2, pair electrons return from 
region III when 7 is as small as 100, because the model re- 
quires that potential variation is of the order of 1/8. Since 
<H/l > 7 (namely j+ > 0) in our model to yield the asym- 
metry of the potential, a larger 7 leads to a larger <f>M- 
As mentioned in Sect. 3.1, </>m should be ~ 1/8. This is 
the reason why we have chosen 7 = 500. When we adopt 
a larger value of 7, the numerical results show that Ef 
becomes smaller. Fig. 4 shows the potential and charge 
densities for one of our results, Case 1. 

Table 1 indicates that a smaller </>m and a larger j p 
lead to a larger amplitude of Ep, as discussed in Sect. 3.1. 
In all our examples the PPF (0 = and I = 0) is in region 
II. The position where the electric field becomes strongest 
is on the boundary of regions I and II. Since the proton 
acceleration is slower than the electron acceleration, the 
boundary is far inside the PPF. For Case 1 the effective 
charge density becomes zero at I — —800 where <f> = —2 x 
10 4 . At this point the amplitude of E becomes maximum. 
In Table 1 we also tabulate the maximum electric field Em 
numerically obtained. The value of Em rather than Ef is 
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Fig. 4. Potential (upper) and charge densities (lower) for 
Case 1. In the lower figure, the thick solid line is the effec- 
tive total charge density. Thin solid line, dashed lines, and 
dash-dotted line are the charge densities due to protons 
(— jp/|/3 p |), injected pairs (/?_ < and p+ > 0), and the 
constant contribution (— (J — jo) = —0.55), respectively. 

the electric field screened in our model. The results show 
that Em is close to the value in Eq. (|25H . 

Region I (I < —800) is not shown in Fig. 4. The ef- 
fective charge density is already positive at the left edge 
(/ = —50) of this figure. The density of positrons initially 
grows outside the PPF, but begins to decrease around 
the point where the positrons injected at the PPF turn 
around. Although the total charge density decreases once 
from this point, the charge density of protons immedi- 
ately makes the total charge density increase again. In 
regions III and IV the densities of both pair-electrons and 
positrons increase while the charge density due to pairs is 
maintaining, p + + p_ ~ The potential of the valley 
for Case 1 is ~ 600 and high enough. 

In our examples cv^ turns out to be 10 -10 -3 , which 
is much smaller than the value required in SMT02. In 
our model the number of electron-positron pairs a beam- 
electron creates per unit time at the PPF is Epa^LOp, 

which becomes ~ 2x ^0 s Tq^B^^ 2 s _1 for Case 1. On 
the other hand, the photon number a primary beam 
electron emits per unit time via curvature radiation is 
2e 2 T/3rh ~ 10 7 r 6 (r/10 7 cm)" 1 s -1 , where r is the ra- 
dius of curvature of the field line. The resultant number 
of pairs created from these photons will be of the same or- 
der. Although we have adopted a very simple model and 
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Table 1. Model parameters and numerical results for E-p and Em 



j P 3+ <I>m j - jo \Ef\ \E m \ 

Case 1 O30 025 1000 055 3.2 x 10~ 4 Isl 25/T 

Case 2 0.30 0.25 600 0.55 1.6 x 10~ 3 21.2 29.3 

Case 3 0.90 0.70 1000 1.60 5.4 x 10" 4 32.9 46.2 

Case 4 0.90 0.70 600 1.60 2.7 x 10" 3 36.8 52.0 



estimate the numbers roughly, these numbers are rather 
close. The required pair creation rate in our examples may 
be consistent with the conventional pair creation models. 

In order to keep the valley shallow enough, our model 
requires that the negative charge due to pairs flows into re- 
gion IV, which leads to the negative charge density of — j+ . 
Since this charge density is very large compared with p; n , 
the total effective charge density remains negative in the 
valley. Therefore, outside the valley the potential grows 
again. Since the charge density due to protons is a func- 
tion of only </>, this charge density dominates again as the 
potential grows, and the potential again begins to decrease 
at (j) = (f>M again. Electron-positron pairs injected after the 
peak flow outwards almost at light speed (the lowest 7_ at 
the valley is ~ 100 for Case 1), so that these pairs do not 
contribute to the charge density. Reiterating this process, 
a space charge density wave continues outside the valley. 
This wave is essentially the same as the wave observed in 
the numerical solutions in SMT98. Of course, the average 
effective charge density is zero in this region. 

4. CONCLUSIONS AND DISCUSSION 

We have shown a possibility that the proton counterflow 
and injected electron-positron pairs screen a strong elec- 
tric field expected in the polar cap model. In our numerical 
examples the required pair creation rate is consistent with 
the prediction of the conventional pair creation models. A 
space charge density wave appears outside the PPF in our 
model. In order to screen a stronger electric field, a larger 
proton current (j p ~ j — j + — jo) and a smaller potential 
difference between the PPF and the peak, 4>m, are favor- 
able. The values of <f>M and the Lorentz factor of pairs at 
injection 7 should be the same order. 

In addition, the amplitude of the charge density wave 
should be small enough to keep pair-electrons from re- 
turning from outside the screening point. The currents of 
returning positrons and protons should be comparable, in 
order to achieve both a small wave amplitude and screen- 
ing of a large electric field. Since a smaller <f>M requires 
a higher pair creation rate, the parameter is limited to 
<t>M ~ O(10 3 ). As for the pair Lorentz factor, 7 ~ 500 
is the best choice. We found that the non-dimensional 
electric field our model can screen is a few times ten at 
most. To attain the Lorentz factor of the primary electron- 
beam, r ~ 10 6 , such a value of the electric field requires 
Pin < 0.01, which means that the screening height is 
greater than 10 4 cm. 

In our model the flux of positrons and protons falling 
onto the star is a few times the GJ flux. The particles 



falling onto the star heat up the polar caps. The X-ray 
thermal luminosity from the neutro n star investigated by 
ROSAT ( Bec ker fc Triimper lll997l) does not indicate that 
there is a significant number of particles falling onto the 
star. Therefore, the fraction of the active region in which 
the electron beam is accelerated in the whole polar cap 
region should be small enough in this model. 

Our model predicts the appearance of a space charge 
density wave after the screening. Particles flow forward 
while accelerating and decelerating in turn. Though the 
spatial average of the effective charge density is zero, the 
local charge density is finite. In more realistic cases, the 
potential wave may not be steady, but oscillate with time. 
The stability of the potential structure may be determined 
by the global dynamics in the magnetosphere, which is be- 
yond the scope of this paper. However, the existence of the 
wave is interesting when we consider the pulsar emission 
mechanism. The wave may induce the charge bunching 
that causes coherent radio emission. 

If most ions in the co-rotating magnetosphere are Fe 
irons, the counterflow could be iron flow. In this case the 
scenario depends on the degree of ionization. If the mass 
and charge of the Fe irons are mp e and qe, respectively, 
5 should be replaced by 8' = q(m p / mp e )8 in our model. 
The smaller q becomes, the larger the typical amplitude 
of the potential variation becomes. Thus we need to adopt 
a larger value of 7 to avoid electrons returning. Even if Fe 
irons are fully ionized (q = 26), 7 is required to be more 
than ~ 1000 in our numerical estimate, while E-p becomes 
larger than in the proton model. 
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